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ABSTRACT 

The central 2x1 kpc of the starburst galaxy NGC 253 has been imaged 
using the Submillimeter Array at a 60 pc resolution in the J— 2-1 transitions of 
12 CO, 13 CO, and C 18 as well as in the 1.3 mm continuum. Molecular gas and 
dust are mainly in the circumnuclear disk of ~500 pc radius, with warm (~40 
K) and high area-filling factor gas in its central part. Two gas shells or cavities 
have been discovered in the circumnuclear disk. They have ~100 pc diameters 
and have large velocity widths of 80-100 km s , suggestive of expansion at ~50 
km s _1 . Modeled as an expanding bubble, each shell has an age of ~0.5 Myr 
and needed kinetic energy of ~1 x 10 46 J as well as mean mechanical luminosity 
of ~1 x 10 33 W for its formation. The large energy allows each to be called 
a superb ubble. A ~10 6 M e super star cluster can provide the luminosity, and 
could be a building block of the nuclear starburst in NGC 253. Alternatively, a 
hypernova can also be the main source of energy for each superbubble, not only 
because it can provide the mechanical energy and luminosity but also because 
the estimated rate of superbubble formation and that of hypernova explosions 
are comparable. Our observations indicate that the circumnuclear molecular 
disk harboring the starburst is highly disturbed on 100 pc or smaller scales, 
presumably by individual young clusters and stellar explosions, in addition to 
globally disturbed in the form of the well-known superwind. 
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1. Introduction 

The center of the nearby spiral galaxy NGC 253 is one of the best places to study 
starburst activity in detail because of its proximity (D = 3.5 Mpc; Rekola et al. 2005). A 
burst of star formation has been taking place in the central ~500 pc for a few 10 7 years (Rieke 
et al. 1980; Engelbracht et al. 1998), producing the far-infrared luminosity of ~3 x 10 10 L Q 
in the shroud of dust (Telesco & Harper 1980). At least part of the star formation has 
occurred in the form of compact star clusters with a few parsec radii and ~10 5 - 10 6 M Q 
masses (Watson et al. 1996; Keto et al. 1999). The young luminous clusters called super 
star clusters have been suspected to be globular clusters in formation. The starburst has 
resulted in the high rate of supernova explosions, ~0.1 yr^ 1 , in the nuclear region (Rieke et 
al. 1988; van Buren & Greenhouse 1994; Ulvestad & Antonucci 1997). The rapid injection of 
energy into the interstellar medium has caused a kiloparsec-scale outflow of hot gas, called 
a superwind, in the direction perpendicular to the galactic disk (Ulrich 1978; Fabbiano and 
Trinchieri 1984; Tomisaka & Ikeuchi 1988; Heckman et al. 1990). 

The fuel for the starburst is the abundant molecular gas in the central kiloparsec, appar- 
ently accumulated with the help of a 7 kpc-long stellar bar (Scoville et al. 1985; Mauersberger 
et al. 1996). Higher-resolution imaging with millimeter arrays showed that molecular gas 
in the starburst region is distributed in a bar-like ~0.7 x 0.2 kpc structure containing sev- 
eral dense clumps, presumably formed in response to the dynamical perturbation from the 
stellar bar (Canzian et al. 1988; Carlstrom 1990; Paglione et al. 1995; Peng et al. 1996; 
Garcia-Burillo et al. 2000; Paglione et al. 2004). The properties of the molecular gas in the 
region are of great importance to further understand the mechanisms and evolution of the 
starburst, because the gas properties affect the formation of stars and are affected by the 
newly formed stars. 

To further study this nearby archetypal starburst, we observed the molecular gas and 
dust in the nucleus of NGC 253 at 1.3 mm using the Submillimeter Array (SMA) 1 on the 
summit of Mauna Kea, Hawaii. This new telescope at 4080 m operates above most of 
the water in the atmosphere, and allowed us to image with high sensitivity and resolution, 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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revealing new features in the starburst region. In this paper, we report the SMA observations 
of the nuclear region with special emphasis on our discovery of two molecular shells that can 
be best modeled as molecular superbubbles. The superbubbles, which are the first ones found 
in NGC 253, are suggested to be the end products of active star formation. We introduce 
our SMA observations in §2, describe the overall properties of gas and dust in the central 
kiloparsec in §3, derive the properties of the expanding molecular superbubbles in §4, discuss 
their formation mechanism and evolution in §5, and summarize our findings in §6. 

We adopt the above-mentioned distance of 3.5 Mpc (1" = 17 pc) for the galaxy, and scale 
the distance-dependent values in the literature accordingly. Also adopted are the following 
parameters; 78° for the inclination of the galaxy (Pence 1980), 51° for the position angle 
(P.A.) of the approaching line of nodes (Burbidge et al. 1962; Jarrett et al. 2003), and ~68° 
for the P.A. of the large scale stellar bar in the galaxy (Scoville et al. 1985). The northwest 
side of the galaxy is the near side to us, assuming that the spiral arms in the galaxy are 
trailing. Figure 1 helps readers to grasp the configuration and to locate the region we discuss 
in this paper. The center of the galaxy is assumed to be at a=00 h 47 m 33 i 517, 5=— 25°17'17'.'l 
(J2000), as has been usually adopted in the literature. The position is determined from the 
brightest compact radio source in the galaxy, possibly an AGN (Turner & Ho 1985; Ulvestad 
& Antonucci 1997). The radio source was suggested to coincide with the galactic center at 
sub-arcsecond accuracy on the basis of the morphology in the infrared and the kinematics 
of dense molecular gas (Kalas & Wynn- Williams 1994; Garcia-Burillo et al. 2000). 



2. Observations and Data Reduction 

The central 2' x 1' (~2 x 1 kpc on the sky) of NGC 253 was observed with the Sub- 
millimeter Array on the 19th and 21st November 2003, during the commissioning of the 
new submillimeter interferometer (see Ho et al. 2004, for description of the telescope). The 
eight 6-meter antennas of the array were in a compact configuration that provided projected 
baselines ranging from 7 m to 75 m, which correspond to spatial frequencies from 4" to 40". 
We observed three positions 26" apart from each other; the interval is the half width at half 
maximum (HWHM) of the primary beam. The central position was at the galactic center 
and the two outer positions were at position angles of 55° and 235°, almost on the major axis 
of the galaxy. The quasar J0132— 169 was observed every 30 minutes for gain calibration, 
Jupiter and Saturn were observed for passband calibration, and Uranus was used for absolute 
flux calibration. Interferometric pointing checks were made toward the quasar J2258— 279 a 
few times before and during the observing run, and toward Uranus before flux calibration. 
Pointing offsets, which typically changed by a few arcsec, were updated about an hour after 
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each pointing measurement because J2258— 279 precedes the galaxy by two hours. The SIS 
receivers were tuned to simultaneously observe three CO lines in the 1 mm band: 12 CO(2-l) 
in the upper sideband, and 13 CO(2-l) and C 18 0(2-l) in the lower sideband. The median 
double-side-band system temperature was 200 K toward the galaxy. The digital correlators 
having a 2 GHz bandwidth were configured for 0.8125 MHz spectral resolution. We lost 
one antenna in the first observing run and lost 12 CO data from two antennas in the second 
observing run because of an antenna and correlator fault, respectively. The total integration 
time on NGC 253 was 7.5 hours. 

Data reduction was made using MIR, which is an IDL version of MMA (Scoville et al. 
1993a), MIRIAD (Sault et al. 1995), and the NRAO MPS package (Bridle & Greisen 1994). 
The data were calibrated for the passband, corrected for atmospheric absorption using the 
system temperatures, and then calibrated for the complex gain using the calibrator data. 
The flux density of the calibrator J0132-169 was estimated to be 1.50 Jy in LSB (220 GHz) 
and 1.44 Jy in USB (230 GHz) from the comparison with Uranus, for which we adopted the 
brightness temperature given by Griffin & Orton (1993). Continuum data were obtained 
from channels that do not have line emission; the resulting data have an effective bandwidth 
of 0.7 GHz in each sideband. The continuum was then subtracted from the visibility data of 
the line emission. Channel maps for the line emission were made with 5 km s -1 resolution 
for 12 CO and 10 km s" 1 resolution for 13 CO and C 18 0. The rms noise is 124 mJy beam -1 
(0.26 K), 44 mJy beam" 1 (0.10 K), and 39 mJy beam" 1 (0.08 K) in the 12 CO, 13 CO, and 
C 18 channel maps, respectively, and is 4 mJy beam -1 in the continuum map made by 
combining the two sidebands. The synthesized beam is almost the same among the lines and 
continuum, ranging from 4'/3 x 2'/5 to 4'/8 x 2'/7. The spatial resolution at the distance of the 
galaxy is about 60 pc. We mosaiced the images using the so-called joint approach, in which 
maps from individual pointing positions were combined first in order to use information 
from all pointing positions in the subsequent deconvolution. Comparison with single-dish 
observations indicates that our interferometric observations recovered about 90 % of the total 
flux in each line as well as in continuum. Specifically, we detected, 91%, 116 %, and 89 % 
of the single-dish integrated intensities of the 12 CO, 13 CO, and C 18 emission, respectively, 
obtained in the central 23" of the galactic center by Harrison et al. (1999). We also detected 
84 % of the single-dish flux density of 1.3 mm continuum that was obtained by Kriigel 
et al. (1990) and is corrected by 30 % for line contamination (Mauersberger et al. 1996). 
Velocities in this paper are with respect to the Local Standard of Rest (LSR) and are defined 
in the radio convention. Unless otherwise noted, maps are shown without correcting for the 
attenuation by the mosaiced primary beam, though all the flux and intensity measurements 
were made after the correction. 
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3. Overall Properties of Gas and Dust in the Central kiloparsec 

The overall distribution of molecular gas and dust in the starburst nucleus of NGC 253 
is shown in the maps of the CO lines and 1.3 mm continuum in Fig. 2. The main component 
is a 50" (0.9 kpc) long linear feature elongated in the position angle of about 58°. The feature 
was first recognized as a molecular gas bar by Canzian et al. (1988) and was subsequently 
modeled as a highly inclined disk near the center of the 7 kpc stellar bar (e.g., Peng et al. 
1996; Garcia-Burillo et al. 2000; Paglione et al. 2004). We simply call the gas structure a 
(circum) nuclear gas disk. The nuclear gas disk extends from about —20" to +30" in the R.A. 
offset in the 12 CO map. Weaker CO emission extends further in the direction of the stellar 
bar. 



3.1. Continuum Emission 

The total flux density of the 1.3 mm continuum is 1.0 ±0.1 Jy in our data, whereas 
that in previous 3.0-3.5 mm observations is 0.25-0.35 Jy (Carlstrom 1990; Peng et al. 
1996; Garcia-Burillo et al. 2000). The 3 mm continuum contains roughly equal amounts of 
synchrotron, free-free, and thermal dust emission (Garcia-Burillo et al. 2000). The increased 
flux density at 1 mm must be from the dust emission because it has a positive spectral slope 
as a function of frequency unlike synchrotron and optically thin free-free emission. Thus 
most (> 80 %) of the 1.3 mm continuum must be thermal emission from dust. 

The continuum emission peaks at a=00 h 47 m 33!18, 5=-25°17'17"4 (J2000) and is elon- 
gated along the nuclear gas disk. The position of the peak agrees, within ~1/10 of our 
synthesized beam, with the adopted galactic center, where the brightest centimeter source 
in the galaxy resides. The radio source called TH2 has been suggested to be an AGN (Turner 
& Ho 1985; Mohan et al. 2002). The millimeter continuum shows a possible secondary peak 
at a=00 h 47 m 32:96, 5=-25°17'19"9 (J2000), which is about 4" southwest of the primary peak 
and is within Of! A from another radio source (called TH8) listed by Turner & Ho (1985). In 
the mid-infrared, the brightest peak in the nuclear region is a super star cluster 3" southwest 
of the galactic center (Keto et al. 1999; Boker et al. 1998). The cluster is located between 
the two millimeter peaks. 

3.2. Molecular Line Emission 

It is known that the bulk of 12 CO(2-l) emission is optically thick while the 13 CO(2-l) 
and C 18 0(2-l) lines are optically thin in the central region of NGC 253, with the mean 
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optical depths in 12 CO(2-l) and 13 CO(2-l) being -3 and ~0.08 respectively (Henkel et 
al. 1993; Harrison et al. 1999). These opacities are based on single-dish observations at 
13"-23" resolutions. The 13 CO line is still optically thin at the higher spatial resolution 
of our data. The maximum value of the 13 CO/ 12 CO brightness temperature ratio is found 
to be 0.4 in our data, after correcting for the slightly different spatial resolution and pri- 
mary beam attenuation. This ratio corresponds to the 13 CO(2-l) optical depth of 0.5, i.e., 
max(r[ 13 CO(2-l)]) ~ 0.5, under the assumptions of LTE and the same beam filling factor for 
the two lines. While 13 CO(2-l) is still optically thin even at the most opaque regions at our 
current resolution, 12 CO(2-l) has a very high opacity in the same regions and suffers from 
self absorption; max(r[ 12 CO(2-l)]) ~ 20, assuming 12 C/ 13 C ~ 40 estimated for NGC 253 by 
Henkel et al. (1993). The same analysis for C 18 results in max(r[C 18 0(2-l)]) ~ 0.15. 

Molecular gas in the nuclear gas disk is warm. For example, the peak brightness tem- 
perature observed in our 12 CO(2-l) data cube is 36 K; the physical temperature that would 
produce the same intensity with the Planck formula is 41 K. The peak intensity is signifi- 
cantly higher than that of the same line in the central kiloparsec of our Galaxy. The peak 
brightness temperature there is 11 K when the 12 CO(2-l) observations by Sawada et al. 
(2001) are convolved to our resolution. The presence of warm gas inferred from the peak 
brightness temperature is consistent with previous observations of NGC 253. The excitation 
analyses of multiple transitions of various molecular lines, such as CO, HCN, and NH 3 , have 
indicated that the bulk of molecular gas is warm, ~50 - 120 K, in the central ~300 pc of 
the galaxy (Glisten et al. 1993; Jackson et al. 1995; Mauersberger et al. 2003; Bradford et 
al. 2003). 



3.3. Structure of the Circumnuclear Gas Disk 

We infer three properties of the circumnuclear gas disk from our data, namely the flaring 
of low density gas from the galactic plane, a high area-filling factor of molecular gas toward 
the center of the disk, and the distinct characteristics of the central 100 pc. 

The flaring of low density gas is inferred from the observation that the 12 CO distribution 
looks wider compared to its length than the 13 CO and C 18 distributions (see Fig. 2). The 
observation was confirmed by fitting an elliptical Gaussian to our data. The deconvolved 
minor-to-major axial ratio is 0.24 in 12 CO, while it is 0.19 and 0.18 for 13 CO and C 18 0, 
respectively. The fitting used maps made by simply summing channel maps rather than the 
moment maps in Fig. 2, in order to avoid the bias caused by the intensity cutoff in the 
moment analysis. The attenuation of the primary beam was corrected before the fitting. 
The fractional uncertainty of the axial ratio is largest for 12 CO and is 2 % for the line, 
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apparently because its distribution is dumbbell-shaped to some extent. For comparison, the 
axial ratio expected from the inclination of the galaxy is 0.20 for a thin uniform disk. The 
broadening in 12 CO may be due to the presence of tenuous molecular gas above and below 
the plane of the galaxy, possibly as part of the superwind. If the flared gas has a density 
slightly below the critical density of the CO(2-l) line (i.e., if uh 2 ~ 10 3 - 10 4 cm" 3 ), it could 
effectively emit only in 12 CO with the help of photon trapping. Entrainment of ambient 
ISM by the superwind has been suggested from extended PAH emission along the minor 
axis (Tacconi-Garman et al. 2005). 

The high area-filling factor of molecular gas in the central part of the nuclear disk is 
implied by the small offset of the 12 CO ridge toward the southeast from the 13 CO and C 18 
ridges as well as from the continuum peak. Figure 3 shows the intensity profiles of those 
lines and continuum through the continuum peak at the position angle of 238°. It confirms 
the offset of / .'5 between 12 CO and its optically thin isotopomers. The offset is unlikely to 
be an artifact of a relative phase-calibration error between the two sidebands, judging from 
the agreement of USB and LSB continuum profiles in the lower panel of Fig. 3. 

A plausible model for the offset is that the molecular disk with a finite thickness has 
a high- area filling factor in its central part. In this model, we see the near-side surface of 
the disk in 12 CO while we see through the disk in optically-thin emission. For a disk of 
half-height h and inclination i, the offset between 12 CO and optically thin tracers will be 
h sin %. The observed offset thus corresponds to the 20 pc thickness (—2h) of the nuclear disk 
measured at its 12 CO photosphere. The disk can be made of gas clumps, but they must be 
so crowded in sky position and in line-of-sight velocity that clump overlap is significant near 
the center of the disk. 

The high area-filling factor (and hence little beam dilution) in the center of the disk 
suggests that the peak brightness temperature of ~40 K is a good estimate for the mean 
physical temperature of molecular gas at the photosphere. For example, it is unlikely that 
there are only small clouds of 200 K filling 1/5 of the beam area. The model also suggests 
that a central hole in the circumnuclear disk, if any, is insignificant at our current resolution. 
If, for example, the nuclear disk were devoid of molecular gas in the central few 100 pc as is 
often observed in the center of barred galaxies (e.g., Sakamoto et al. 2004), then that should 
show as a clear central depression in optically thin emission. Without such a feature, the 
nuclear disk must have molecular gas (with a high-area filling factor) in its center, though 
the gas disk could still have a ring-like region of high density around the galactic center as 
suggested from some observations (Israel et al. 1995; Garcia-Burillo et al. 2000; Weaver et 
al. 2002). 

The distinct characteristics of the central 100 pc of the nuclear gas disk is indicated by 



- 8- 



the markedly different spatial distribution of the 1.3 mm continuum and that of 13 C0 and 
C 18 integrated intensities. The former is more strongly peaked toward the galactic center 
than the latter. This may be because free-free emission, which is less than 20 % of the 1.3 
mm continuum, is strongly peaked toward the galactic center. Indeed the 1.3 cm radio map 
of Ulvestad & Antonucci (1997, see their Fig. 2) would look like the 1.3 mm map near the 
galactic center if convolved to the same resolution; the brightest peak is at TH2 and the main 
extension is to the southwest. If the free-free emission is the cause of the strong continuum 
peak at the galactic center, then that means the star formation is most intense there. Another 
class of possible reasons for the discrepancy between continuum and 13 CO as well as C 18 is 
the variation of the properties of the neutral ISM near the galactic center. The intensity of 
1.3 mm dust emission should approximately be the product of dust mass opacity coefficient, 
the dust column density, and the temperature of the dust, or /dust, 1.3 mm «i.3 mm Sd us t Tkust- 
The intensity of an optically thin CO (2-1) line is approximately proportional to the column 
density of molecular gas and the fractional abundance of the CO molecule with respect to 
H 2 , or ico(2-i) oc Sh 2 ([C0]/[H 2 ]), if the density and temperature of the gas are in the range 
of 40-200 K and 10 3 - 10 5 5 cm" 3 , respectively (Wild et al. 1992). The bulk of molecular gas 
in the nuclear disk satisfies the condition (Harrison et al. 1999; Paglione et al. 2004). Thus, 
if non-thermal emission is not strongly peaked toward the galactic center, possible reasons 
for the higher continuum-to- 13 CO (or C 18 0) ratio within the central 100 pc of the nuclear 
gas disk include a higher temperature of the ISM, higher dust-to-gas ratio, and the variation 
of the dust opacity coefficient or that of the CO isotopomer abundance. Each of these cases 
again suggests that the central 100 pc has distinct characteristics in the nuclear gas disk, 
likely due to the intense star formation there. 



4. Expanding Molecular Superbubbles 

The 12 CO channel maps in Fig. 4 show two shell-like features in the nuclear gas disk. A 
similar shell-like feature of comparable size has been observed in the starburst nucleus of M82, 
and has been interpreted as an expanding molecular superb ubble (Neininger et al. 1998; Wills 
et al. 1999; Weifi et al. 1999; Matsushita et al. 2000; Weifi et al. 2001; Matsushita et al. 2005). 
We hence call the one around Vlsr ~ 380 km s _1 SB1 and another one around Vlsr ~ 150 
km s _1 SB2, under the working assumption that they are molecular superbubbles. Both 
features are also seen in our 13 CO data, and SB1 is barely visible even in the C 18 data. 
SB1 is at a kink of the almost linear gas distribution (see Fig. 2), and may be responsible for 
it as we discuss below. The kink was already seen in the first millimeter-interferometer image 
of the nuclear gas disk (Canzian et al. 1988), though the shell feature there and another one 
on the other side of the nucleus were not recognizable in previous observations because of the 
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lack of sensitivity or spatial resolution or both. In the following, we derive the parameters 
of the shells from our data and summarize them in Table 1. 



4.1. Size, Shape, and Location 

The size of the shells, estimated from the channel maps and shown in the maps, is about 
160 x 100 pc. The shells look elongated in the southeast-northwest direction, roughly along 
the minor axis of the galaxy. The exact size and shape of the features are, however, difficult 
to determine at the current resolution. Part of the elongation must be due to the elongation 
of the synthesized beam in the same direction. Also, it is almost certain that the shells 
have more complex shapes than an ellipse in the plane of the sky. For example, the shells 
may not be closed but have openings near its major axis as hinted in the channel maps, 
such as 350, 375, and 410 km s _1 maps of SB1. These openings might indicate a conical 
or cylindrical shape elongated along the minor axis of the galaxy, as would be expected for 
galactic bubbles. It is also noted that the channels from around 370 km s _1 to 400 km s _1 
show an emission tongue that appears to emanate from the area of SB1 by up to 10" (170 
pc) in the minor axis direction. Given the large inclination of NGC 253, the tongue emission 
is probably extended in the direction perpendicular to the disk rather than within the plane 
of the galaxy. 

The centroid of each shell was determined from the local minimum in the map that was 
made by integrating over the central velocities of each shell (Fig. 5). Their distances on 
the sky from the center of the galaxy are 11" and 22" for SB1 and SB2, respectively. Their 
deprojected distances from the galactic center would be 190 pc for SB1 and 670 pc for SB2 
if the centroids are on the midplane of the galaxy. SB2 could be at a smaller galactocentric 
radius if its center is above (i.e., in front of) the midplane of the galaxy; a 50 pc vertical offset 
would place the centroid at the galactocentric radius of 500 pc. This may be a more plausible 
configuration because otherwise the circumnuclear disk would be lopsidedly elongated in the 
direction of SB2. Note, however, that the geometrical center of an ISM bubble does not 
necessarily coincide with the location of its driving source, because a bubble expands faster 
in the direction of lower ambient pressure. For SB1, any vertical offset would make its 
galactocentric radius larger than 190 pc, because the center of the feature is almost on the 
line of nodes. 
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4.2. Kinematics, Ages, and Orbits 

The full- velocity widths in which the shells are visible in the channel maps are about 100 
and 80 km s _1 for SB1 and SB2, respectively. The velocity structure of SB1 is clearly seen 
as an elliptical hole or a depression in the 12 CO and 13 CO position- velocity (p-v) diagrams 
in Fig. 6. The hole structure is similar to the one seen in the p-v diagrams of M82 at the 
location of the molecular superbubble in the galaxy (see Fig. 2 of WeiB et al. 1999). The 
shell feature SB2 is less evident in p-v diagrams; a partial hole is marginally visible in 13 CO 
though it does not show up as a hole in the 12 CO p-v diagram (Fig. 7). A hole in a p-v 
diagram suggests that the shell is around a three-dimensional cavity that is either expanding 
or contracting. The latter is highly unlikely because a coherently contracting shell of ~100 
pc size is difficult to create, whereas an expanding shell can be produced around a point 
explosion. The expansion caused by an explosion would be faster along the minor axis of the 
galaxy because of the decreasing pressure away from the disk. The elongation of the bubbles 
in Figs. 4 and 5, if confirmed, would be consistent with this expectation. 

Assuming that the shells are expanding bubbles, the kinematical age of each shell is 
estimated to be ps 0.5 Myr. This is derived from tkin ~ cxR/v, where R is the radius and 
v is the expansion velocity. We use the semi-minor axis and half velocity width of each 
bubble for R and v respectively, so that both are the ones in the plane of the nearly edge-on 
galaxy. Cylindrical symmetry of each feature is assumed. The parameter a is to account 
for non-linear expansion in the form of R oc t a . For example, a is 1 for free expansion, 3/5 
for adiabatic expansion due to constant energy injection, 2/5 for non-radiative expansion 
due to a single explosion (the Sedov- Taylor phase), 2/7 for the pressure-driven expansion 
of a shell around a hot bubble, and 1/4 for momentum-conserved snowplow expansion (see, 
e.g., Ostriker & McKee 1988); all the explosions are assumed to take place in a uniform 
interstellar medium with negligible pressure. The power a could be larger than 3/5 if the 
energy injection is continuous and has been increasing with time, as one might expect from a 
gradually forming star cluster. We adopt a 0.5 as a value in the middle of possible values 
and expect a factor of 2 uncertainty in the estimated age. 

We can also infer orbital properties of the bubbles from their mean velocities and the 
p-v diagrams. There are basically two types of oval orbits for gas and stars in a galaxy 
bar. Namely, the dominant X\ orbits are elongated along the bar, and x 2 orbits elongated 
perpendicular to the bar exist near the galactic center under certain circumstances (c.f., 
Binney & Tremaine 1987). Gas can take intermediate orbits near the transition region 
between the two types of orbits. In the configuration of NGC 253, a gas cloud in an x 2 orbit 
should have larger line-of-sight velocity than that in an x\ orbit, as modeled by Peng et al. 
(1996). The presence of two types of orbits (i.e., velocity components) in the central ~50", 
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where there is the nuclear gas disk, is clearly seen in the single-dish CO p-v diagram of Sorai 
et al. (2000, their Fig. 17(a)). The high-velocity component there is most likely in the x<i 
orbits or in the orbits that are intermediate between x\ and rr 2 , whereas the low- velocity 
component within about 50 km s _1 from the systemic velocity must be in the X\ orbits. The 
bubble SB1 and its progenitor should belong to an X2 orbit because the mean velocity of 
the bubble is close to the terminal velocity in the p-v diagram and about 130 km s _1 away 
from the systemic velocity. The association of SB2 to another x 2 orbit is likely but less 
certain, because the offset from the systemic velocity is 90 km s -1 . The ages of the bubbles 
are several ten times shorter than their orbital periods around the galactic center. Thus, 
differential rotation could not have significantly affected the bubbles. 



The total energy needed to create each bubble is estimated in two ways to be ~10 46 J 
from the size and expansion velocity of the bubble as well as the mean density of the ambient 
gas. First, the numerical simulation of a point explosion in the uniform interstellar medium 
by Chevalier (1974) gives the following formula for the total energy of the single explosion 
that caused the bubble: 



where no is the hydrogen number density of the pre-bubble medium. Alternatively, one 
can use the similarity solution of a bubble caused by continuous energy injection into the 
interstellar medium. The total injected energy, 20% of which goes to the kinetic energy in 
the solution (Weaver et al. 1977), is estimated to be 



where p = lAn mu is the initial mass density of the ISM. The gas density in both formulae 
is the one prior to the bubble formation and is assumed to be uniform. One therefore has 
to estimate the pre-bubble gas density averaged over the volume that the bubble currently 
occupies. We adopt no ~ 70 cm~ 3 and ~ 20 cm -3 for SB1 and SB2, respectively, on the basis 
of the gas column densities estimated for the nuclear gas disk (Wall et al. 1991; Mauersberger 
et al. 1996; Harrison et al. 1999; Sorai et al. 2000). Specifically, we used 12 CO(2-l) integrated 
intensities of 2 x 10 3 and 0.5 x 10 3 K km s _1 for the bubbles, the inclination of the disk of 
78°, NuJIco = 0.3 x 10 20 cm" 2 (K km s -1 ) -1 , and a 100 pc height for the bubbles. The 
kinetic energy needed to create the R ~ 50 pc bubbles is estimated to be 2-3 xlO 46 J for 
SB1 and 0.4-0.6 xlO 46 J for SB2. The smaller numbers are from the similarity solution, in 
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which radiative loss of energy is ignored. The mass of each shell is M she ii ~ 1.3 x 10 6 M 
and 0.4 x 10 6 M for SB1 and SB2, respectively. 

It is noted that the error of the derived energy may be as large as an order of magnitude 
in the worst case. This is because the energy depends almost linearly on the ambient gas 
density, whose uncertainty is no less than a factor of a few, and also depends on many 
parameters that are not in the models, such as non- uniformity of the ambient gas. The 
pressure of the ambient gas is among those ignored in the models even though the pressure 
300 pc from the center of NGC 253 has been estimated to be ~10~ 10 Pa (i.e., P/k ~ 10 7 cm -3 
K; Heckman et al. 1990; Carral et al. 1994), which is about 2 orders of magnitude higher 
than the pressure in the center of our Galaxy (Sawada et al. 2001). The continuous energy- 
injection model gives the internal pressure of a bubble as P in = -^Pov 2 , and each bubble has 
~5 times higher internal pressure than external pressure according to the formula. Thus the 
bubbles probably have not been affected significantly by the external pressure yet. However, 
now that there are morphological hints of blowouts and the size of the bubbles is indeed 
comparable to the thickness of the molecular gas layer in galactic centers (Scoville et al. 
1993b, and references therein), the bubbles are likely to be approaching the end of the 
pressure-driven expansion in the disk plane of the galaxy. 

The energy needed to create each bubble, ~10 46 J, corresponds to the kinetic energy 
released by ~100 supernovae (SNe), or ~1— 10 hypernovae (HNe) each of which releases 10- 
100 times more kinetic energy than a supernova (Iwamoto et al. 1998; Nomoto et al. 2004). 
This suggests that, unless the energy came mostly from a hypernova, each of the expanding 
bubbles was created by successive or continuous injection of energy rather than by a single 
explosion. The formation of a large ISM bubble by continuous energy injection is consistent 
with the model of superbubbles, which are created by stellar winds and supernova explosions 
of massive stars in a star cluster (Tomisaka et al. 1981; McCray & Kafatos 1987). It is also 
reasonable to call an ISM bubble caused mainly from a hypernova explosion a superbubble, 
because the involved energy matches that of multiple supernovae in the original superbubble 
model. Thus, SB1 and SB2 can be called superbubbles in terms of energetics. 

4.4. Association with Features in other Wavelengths 

There is only one compact radio source inside the bubbles according to the list of 
centimeter- wave continuum sources of Ulvestad & Antonucci (1997), even though about 60 
compact sources most likely comprised of supernova remnants and HII regions have been 
detected in the nuclear region. Most of the radio sources are concentrated around the 
starburst nucleus between the two bubbles. The radio source listed as 5.17—45.4 is within 
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1" from the center of SB1 (see Fig. 5) and has a 6 cm flux density of 0.32 mJy (Antonucci & 
Ulvestad 1988). The source by no means stands out among the radio sources in the nuclear 
region, nor is it very prominent in terms of its absolute radio power. The radio power of the 
source, P^m = 5 x 10 17 W Hz -1 , is equivalent to that of Cas A, a single supernova remnant 
in the Galaxy. The lack of a significant population of centimeter sources in the bubbles, 
however, does not mean that there have been few supernovae. This is because the VLA 
observations can detect a supernova remnant for only several 100 years, as can be estimated 
from the supernova rate, 0.1 yr -1 , and the number of radio sources detected in the starburst 
region. For a comparison, the mean interval between supernovae in the bubble would be 
5000 yr to have 100 SNe in 0.5 Myr. The bubbles do not coincide with known OH and H 2 
masers, which are concentrated to the central 200 pc of the starburst nucleus (Frayer et al. 
1998; Henkel et al. 2004). 

There is a faint spot near the center of SB1 in the HST image of the [Fell] 1.644 /im 
emission line (Alonso-Herrero et al. 2003, Fig. 5b). The spot is one of the compact sources 
that show spatial correspondence with a centimeter radio source. Compact [Fell] sources 
in a starburst are thought to trace supernova remnants for ~10 4 yr. In the HST study, 
the median [Fell] luminosity of the compact sources is two orders of magnitude larger than 
that of galactic supernova remnants. Thus the sources may well be clusters having multiple 
supernova remnants, though higher density of the ISM could also boost the line luminosity 
(Morel et al. 2002). In X-rays, inspection of Chandra images reported in Strickland et al. 
(2000) and Weaver et al. (2002) did not reveal any obvious source in the bubbles. 

The two superbubbles are around the waist of the superwind, which is a kiloparsec-scale 
conical outflow of ionized gas from the starburst region into the direction perpendicular to 
the galactic disk. The base of the superwind has a diameter of ~20" in X-ray, near-IR H 2 line, 
and synchrotron emission (Strickland et al. 2000; Sugai et al. 2003; Ulvestad & Antonucci 
1997). The bubble SB1, which is 11" from the galactic center, is almost at the surface of 
the biconical outflow. The CO tongue extending northwest from around SB1 corresponds to 
the H 2 spur B of Sugai et al. (2003) and also coincides with the foot of the most prominent 
spur in the 20 cm image of Ulvestad & Antonucci (1997). There are two possible causes for 
the CO tongue if it is a ~200 pc out-of-plane structure as we inferred. One possibility is a 
blowout from the superbubble SB1 accidentally overlapping with the superwind cone, and 
the other is molecular gas entrained by the large-scale superwind. 
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4.5. Gas Properties in and around the Bubbles 

A superbubble is expected to affect the state of the ISM in and around the expanding 
shell. In an attempt to see the effect, we made a position-velocity diagram through SB1 for 
the 13 CO-to- 12 CO temperature ratio. The p-v diagram shown in Fig. 8 indicates that the 
line ratio in the bubble region is lower by a factor of a few than the ratio at the ambient 
non-bubble locations. 

There are a number of possible reasons for the lower ratio at the bubble location. 
The lower ratio can be simply due to the lower column density per unit velocity in the 
bubble, because that reduces 12 CO opacity. Also, the line ratio can be due to the excitation 
conditions of the CO lines. In the bubble, the number density of molecular gas may be below 
the critical density to excite CO to the J=2 level. It would make the excitation of 13 CO 
much harder than that of 12 CO because photon trapping does not help the optically thin 
13 CO. Alternatively, the gas temperature in the bubble can be high enough to move most 
population above J=2, making both 12 CO and 13 CO optically thin, and reducing the intensity 
ratio to the abundance ratio of the two molecules. In addition, selective photodissociation of 
13 CO molecules (van Dishoeck & Black 1988) may result from intense UV radiation in the 
bubble, again reducing the line ratio. Observations of other transitions and molecules would 
help to distinguish between these possibilities. 

5. Evolution of the Superbubbles, and Their Progenitors 

We have seen that the ~100 pc shell features found in the molecular gas around the 
nuclear starburst of NGC 253 can be modeled as expanding superbubbles of ages ~0.5 Myr 
and injected kinetic energies of ~10 46 J. The evolution and progenitors of such superbubbles 
can be inferred from the current data. 



5.1. Young Bubbles and their Evolution 

The kinematical ages of the superbubbles, ~0.5 Myr, are rather short compared to the 
age of the starburst, which is estimated to be 20-30 Myr (Engelbracht et al. 1998). Why are 
older superbubbles missing? 

One of the possible solutions for this mystery is that older bubbles are difficult to find 
because they must have blown out of the galactic plane and have ceased expansion. If a 
superbubble similar to the ones we found has reached a point of blowout, as at least SB1 
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seems to have, then its expansion within the galactic plane loses the thrust from internal 
pressure and becomes a momentum- driven snowplow. The expansion in the plane of the 
galaxy halts (i.e., the expansion velocity reduces by a factor of 5 and becomes comparable to 
the turbulent velocity of the ambient molecular gas) in about 2 Myr. The size of the bubble 
increases only by about a factor of 2 in that time. After that the bubble will no longer have 
a quasi-spherical shape nor significant expansion velocity when viewed from edge-on. Thus 
such a superbubble would be difficult to recognize despite its large size. An older bubble 
would get less recognizable also because its outer wall becomes fragmented and because its 
shape in the galactic plane is deformed from the original circular one owing to the differential 
rotation and non-circular gas motions in the region. If a typical molecular superbubble could 
be detected for only a few Myr with one's resolution and sensitivity, then it may well happen 
that one finds two half-Myr old superbubbles without seeing older ones. 

It is also possible that star formation in the radii where we found the superbubbles has 
been retarded compared to the nuclear starburst inside, reaching the stage of superbubble 
formation only recently. The delay could be due either to the inside-out progression of 
star formation or to the time lag statistically expected from the low superbubble formation 
probability in the outer region. In contrast, the nuclear region with a ~10 7 -yr-old starburst 
may be already too disturbed or have too high ambient pressure to form a new distinct 
superbubble. The energies from stellar winds and supernovae there may be channeled almost 
directly to the superwind. 

We note that the molecular superbubble found in M82 is also very young, having a ~1 
Myr age (WeiB et al. 1999; Matsushita et al. 2000). Being young (only about 1 Myr old) is a 
common character of the molecular superbubbles found so far in two of the nearest starburst 
galaxies. This may well be due to the short time window of detection for superbubbles in 
edge-on galaxies, because both starburst galaxies are seen nearly edge-on. 

5.2. Super Star Clusters 

The large energy and the young age of each superbubble provide a strong constraint on 
the progenitor, because a large mechanical luminosity is required to create the bubble. The 
progenitor of each bubble must have had a mean mechanical luminosity of 1 x 10 33 W in 
the last 0.5 Myr in a region smaller than 100 pc. This luminosity can be provided by a 10 6 
Mq star cluster that has the Salpeter initial mass function, an upper mass cutoff of 100 M Q , 
and solar metallicity, according to the Starburst99 model of Leitherer et al. (1999). In the 
case of instantaneous star formation in the cluster, the luminosity stays roughly constant for 
about 30 Myr, first coming from stellar winds and then from supernovae. The duration of 
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successive supernova explosions, ~30 Myr, is set by the range of the lifetimes of the massive 
stars (~8-lOOM ) that eventually explode as core-collapse supernovae. The model cluster 
at the age of 1 Myr has 5000 O stars, absolute K-band magnitude of Mr ps — 16 mag, and 
bolo metric luminosity of 1 x 1O 9 L , or 3 % of the total far-IR luminosity of the starburst 
and that of the galaxy (Lfir = (3.0-3.5) x 1O 1O L , Telesco & Harper 1980; Sanders et al. 
2003). A star cluster of this magnitude, if it exists, must be a significant building block of 
the starburst in NGC 253. 

The mass of the cluster postulated here agrees with that of the super star clusters found 
in many starburst galaxies including NGC 253 (Whitmore 2003, and references therein), 
though the cluster mass, as well as the numbers derived from it, has uncertainty inherited 
from the uncertainties in the expansion energy and age of the superbubbles. Super star 
clusters are therefore possible driving sources of the superbubbles. The measurements of 
the size and stellar population of the progenitor clusters, if any, are needed to conclude the 
identification. 

At the distance of NGC 253, the cluster's magnitude will be m K 12 mag without 
extinction, according to the Starburst99 model. Extinction toward the center of each bubble 
is estimated from average H 2 column densities to be A K < 5 and < 2 mag for SB1 and SB2 
respectively, though small scale structure in the molecular disk may considerably change 
local extinction in our lines of sight. Thus the clusters would nominally be at 14-17 mag in 
the K band. The bubble locations have i^-band surface brightness of 14-15 mag arcsec -2 
in the 2MASS Large Galaxy Atlas (Jarrett et al. 2003). The 2"-3" resolution of the data, 
however, is insufficient to identify compact sources in the region. The high-resolution near- 
IR images of Alonso-Herrero et al. (2003) do not show a conspicuous source at the bubble 
locations in continuum, though a source is seen in SB1 in the [Fe II] line as we already 
noted. The datasets we found in the literature do not allow us to set a constraint in the K 
band. In 1 mm continuum, the uncertainty in dust configuration in our lines of sight to the 
progenitors makes it difficult to set a constraint. On one hand, a cluster of 1 x 1O 9 L should 
have been detected at ~30 mJy, which is ~7<r in our continuum map, if dust reprocesses the 
cluster luminosity by the same degree as it does to the bulk of the starburst luminosity in 
NGC 253. On the other hand, the conversion to dust emission may be less effective for the 
progenitors considering that they should be in cavities of interstellar medium. 

If the two superbubbles are indeed made by star clusters, then the bubbles must have 
spent only a few % of their energy injection period until now. In other words, the total 
kinetic energy eventually released inside each superbubble through SNe would be ~10 48 J, 
or several 10 times larger than the energy deposited so far. Because a superbubble blows 
out of the disk early, most of the energy will be channeled to the galactic halo as depicted 
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in the galactic chimney model by Norman & Ikeuchi (1989). In a similar starburst in M82, 
several chimneys are indeed seen in deep radio and optical images (Wills et al. 1999). The 
detection of such chimneys in NGC 253, however, would be more difficult because NGC 253 
is not completely edge-on. 

5.3. Hypernovae 

A hypernova can, as we saw, provide most of the mechanical energy needed to form 
each of the superbubbles. It also satisfies the constraints on the mean mechanical luminosity 
and compactness because it is an instantaneous point explosion. Furthermore, it is possible 
from a statistical viewpoint that the observed molecular superbubbles are due to a couple 
of hypernovae. The supernova rate in the nuclear starburst in NGC 253 is estimated to 
be about 0.1 yr _1 (Rieke et al. 1988; van Buren & Greenhouse 1994; Ulvestad & Antonucci 
1997). The hypernova-to-supernova ratio, on the other hand, has been estimated to be 1CT 3 - 
1CT 5 (Paczyhski 1998; Podsiadlowski et al. 2004). Thus the hypernova rate in the starburst 
region of NGC 253 is expected to be in the range of ~10~ 5±1 yr _1 . This is comparable 
to or larger than the superbubble formation rate of 10~ 6 yr^ 1 estimated from the current 
number of superbubbles, namely 2, and the time scale they are recognizable, ~2 Myr. It is 
therefore possible, not only energetically but also from the formation rate argument, that 
the superbubbles were caused by hypernova explosions. Note that there must be 1-100 
hypernova remnants of less than 1 Myr old in the starburst region, regardless of whether the 
two superbubbles were due to hypernovae. 

The total kinetic energy released from the hypernova population is much less than that 
from supernovae, when integrated over the starburst region and over the entire duration of 
the starburst. It is because the larger number of supernovae more than compensates for 
the smaller energy of each explosion. However, supernovae can create a superbubble only 
if they are concentrated both in space and time, as in a star cluster. This constraint does 
not apply to hypernovae so long as the superbubbles that could be created with a single 
hypernova explosion are concerned. Therefore hypernovae can be more significant than 
expected from the total-energy argument as a mechanism to create molecular superbubbles in 
an environment without massive clusters. On the other hand, a hypernova in an environment 
with many massive clusters is more likely to occur in one of those clusters and to cooperate 
with stellar winds and supernovae to form a superbubble around the cluster. 

Observationally, the largest difference between a super star cluster and a hypernova as 
the progenitor of a molecular superbubble is that the former must be still present as a cluster 
of thousands of OB stars. If such a cluster is not found in a superbubble in any wavelengths, 
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then the progenitor is more likely to be a hypernova. 

Regarding the effect on the ISM, a hypernova will have less lasting impact than a super 
star cluster, even though each forms a similar superbubble when compared at the age of ~1 
Myr. The difference is because, as we noted, massive stars in a super star cluster last for a 
few 10 Myr and release ~100 times more kinetic energy than a hypernova. 

6. Discussion 

Our observations revealed many gas properties and features that are most likely caused 
by the starburst. They include the high gas temperature, superbubbles, and molecular gas 
likely entrained by the superwind. 

Regarding the bubbles, we saw that a number of properties of the two shell features 
are consistent with them being expanding superbubbles. To be precise, however, that does 
not prove them to be superbubbles. Cold gas in a barred galaxy has non-circular motion 
and non-axisymmetric distribution due to the bar, and it could show a bubble-like feature 
in a p-v diagram, as pointed out by Wills et al. (2000). We therefore tested if we see the 
shell-like features in gas dynamical models in various barred potentials, by taking the models 
of Athanassoula (1992) and plotting them in the same way as NGC 253 was observed. The 
bubble features in the sky-plane maps and p-v diagrams could not be reproduced in the 
models we tested. In particular, the bar models have an inherent difficulty to explain two 
shell-like features whose positions are not symmetric about the galactic center, because bar- 
induced features tend to be bisymmetric. Some additional factors, such asm ^ 2 modes 
in the bar potential and self-gravity of gas, are required to break the bisymmetry. We have 
not explored such models, because the parameter space becomes vastly larger for them. In 
this regard, there remains an untested possibility that factors unexplored by us conspired 
in NGC 253 to create gas features masquerading as expanding bubbles. We just note that 
such gas features, if confirmed, would be rarer findings than the molecular superbubbles that 
have been reported in several galaxies (see e.g., Handa et al. 1992; Irwin & Sofue 1996; Rand 
2000; Walter et al. 2004, and those already cited for the molecular superbubble in M82). We 
conclude that the current data most favor the superbubble model. Further observations, in 
particular the ones to look for the progenitors or their remnants in various wavelengths, are 
desirable. 

Unusual features in the central molecular disk in NGC 253 are not limited to the couple 
of gas bubbles. The 12 CO position-velocity diagram shown in Fig. 6 is highly disturbed in 
comparison with those in the collection of CO p-v diagrams in Sakamoto et al. (1999). For 
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example, the trough at the offset of —11" and the emission components on both of its sides 
having nearly 150 km s^ 1 velocity widths are quite unusual. The gas dynamical simulations 
we mentioned above did not produce the trough feature either. These features may be related 
to the superwind cone around that location or may be due to another superbubble, a very 
young one that would be characterized by its compact size and large velocity width. 

The numerous signs of disturbances in the central disk once again prove NGC 253 to 
be one of the best targets to study starburst. The features we found are of smaller scales 
(i.e., 100 pc or less) than the ones previously known, thanks to the higher resolution of 
our observations. These kinematical and thermal inputs to the local interstellar medium 
should almost certainly influence the subsequent star formation in the region, through the 
formation and destruction of star forming gas clouds. These inputs should therefore affect 
the evolution of the starburst. Details of this process remain to be seen in further studies. 

7. Summary 

We observed molecular gas and dust in the central ~2 kpc of the archetypal starburst 
galaxy NGC 253 at resolutions of ~60 pc and 5-10 km s -1 . Our main conclusions are the 
following: 

1. Molecular gas in the central region is mostly distributed in the nuclear gas disk of ~500 
pc radius. The gas in the disk is warm, with a peak brightness temperature of ~40 K 
in 12 CO(2-l). The area-filling factor of the gas appears high in the center of the disk, 
suggesting the observed brightness temperature to be close to the gas temperature. 
There are hints that low density molecular gas is flaring above the midplane possibly 
due to the superwind, and that the ISM in the central 100 pc has distinctive physical 
conditions. 

2. We have found two molecular shells or cavities of ~100 pc size in the nuclear gas disk. 
These features show velocity structure indicative of expansion at a velocity of ~50 
km s -1 . 

3. We modeled these features as molecular superbubbles. In the model, each feature has 
an age of ~0.5 Myr, involves ~10 6 M & of interstellar medium, and needed the kinetic 
energy of ~1 x 10 46 J for its formation. 

4. The large amount of energy deposited within the short timescale calls for a progenitor 
with a high rate of energy injection, 1 x 10 33 W averaged over the age of the superbub- 
bles. A ~10 6 M Q star cluster of ~1 Myr age would have enough mechanical luminosity 
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to create such a superbubble with the stellar wind and from supernova explosions of 
thousands of OB stars. If such a cluster had a small size of < 10 pc, it would be a 
super star cluster. It could contribute a few percent of the total bolometric luminosity 
of the starburst in NGC 253, and would be a significant building block of the starburst. 

5. Alternatively, a hypernova, which releases 10-100 times more kinetic energy than a 
supernova, could be a significant source of energy for each superbubble. The hypernova 
rate in the starburst region, inferred from the supernova rate there, is of the same order 
of or larger than the formation rate of the superbubbles. Thus a hypernova is a tenable 
progenitor of each superbubble. 

6. Our data suggest that the nuclear starburst in NGC 253 had strong thermal and 
kinematical impacts on the central molecular disk hosting the starburst. In particular, 
our high-resolution observations have revealed various signs of disturbance at 100 pc or 
smaller scales in the gas disk. These may well influence the formation and destruction 
of star forming gas clouds, and thereby affect the evolution of the starburst. 
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Table 1. Molecular Superb ubbles in NGC 253 



parameter 


SB1 


SB2 




(J0047325-251724) 


(J0047346-251709) 


a J2000 


00 h 47 m 32!53 


00 h 47 m 34!64 


<5j2000 


-25 o 17'24'.'0 


-25°17'09'. / 5 


V LS r [km s" 1 ] 


375 


150 


size [arcsec] a 


9x6 


10 x 6 


major axis P. A. a 


135° 


135° 


AV [km s- 1 ] 


100 


80 


4m [Myr] b 


0.5 


0.6 


Rg.c. [pc] c 


190 


< 500 (or 670) 


E [10 46 J] 


2 


0.4 


M shell [1O 6 M ] 


1.3 


0.4 



Approximate size and position angle. See text for details. I" 
= 17 pc. 

b Kinematic age. 

c Distance from the galactic center. See text for the discussion 
on the two possible assumptions about SB2. 
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Fig. 1.— (top) The J - K s color index image of NGC 253 made from the 2MASS data of 
Jarrett et al. (2003). The gray scale is from 0.9 mag to 2.3 mag with darker regions being 
redder, or showing higher extinction by interstellar dust. The dashed ellipse, only a part 
of which is in the frame, is the galaxy's outline at K s = 20 mag. The dotted oval in the 
center shows the field of view of our mosaic CO observations. North is up and east is to the 
left. The northwest side, i.e., the top-right side, is the near side of the galaxy, (bottom) An 
illustration showing the main features in the galaxy. The new dashed ellipse is the stellar 
bar that is 7 kpc long on the plane of the galaxy. The gray thick curves delineate dust lanes 
in the leading edge of the stellar bar as well as the spiral arms starting at the ends of the 
bar. The gray ellipse in the center represents the circumnuclear disk of molecular gas where 
a burst of star formation has been taking place. The dash-dotted line is the line of nodes. 
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Fig. 2.— The starburst nucleus of NGC 253 in 12 CO(2-l), 13 CO(2-l), C 18 0(2-l), and 
1.3 mm continuum. Positional offsets are measured from a=00 h 47 m 33!17, 5=— 25°17 / 17'. / 1 
(J2000), which is the assumed position of the galactic center and is shown as a cross. The 
brightest centimeter-wave source in the galaxy, called TH2, is at this position. The unit of 
intensity is Jy beam -1 km s _1 for the CO integrated intensity maps and Jy beam -1 for the 
continuum map, which is made by averaging the USB and LSB maps. The FWHM of the 
synthesized beam is shown at the bottom left corner of each map, and the mosaiced primary 
beam is indicated by a dotted line at 40 % of its peak. The positions of the molecular 
superbubbles are shown with asterisks, and their approximate sizes are shown with ellipses 
in the 12 CO image. The superbubble called SB1 is to the southwest and the one called SB2 
is to the northeast of the nucleus. 
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Fig. 3.— The intensity profiles of 12 CO(2-l), 13 CO(2-l), C 18 0(2-l), and 1.3 mm continuum 
across the continuum peak along the position angle of —32°, which is approximately the minor 
axis of the nuclear gas disk. The profiles are measured from the maps that are corrected for 
the primary beam attenuation and convolved to the same spatial resolution (the resolution 
of the C 18 map). The CO line maps are made by summing up channel maps rather than by 
moment analysis. The amplitude of each emission is normalized so that the peak intensity is 
unity. Note that each emission is only partially resolved in this direction. The peak signal- 
to-noise ratio is 330, 75, 29, 29 and 39 for 12 CO, 13 CO, C 18 0, and USB and LSB continuum, 
respectively. The 12 CO(2-l) profile is 075 offset from that of 13 CO(2-l) and C 18 0(2-l). 
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Fig. 4. — 12 CO(2-l) channel maps around the superbubble SB1 (top) and SB2 (bottom). The 
n-th contour is at n 1 ' 2 x 0.5 Jy beam -1 (=4<r) for SB1 and n 1,3 x 0.37 Jy beam -1 (=3<r) for SB2. 
The conversion to excess brightness temperature is 1 Jy beam -1 = 2.1 K. The centroid position and 
approximate size of each bubble are marked with a red cross and a red ellipse, respectively. Velocity 
is at the top-right corner of each panel, and the synthesized beam is shown at the bottom-left corner 
of the first panel. 
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Fig. 5.— The molecular superbubbles in NGC 253 in 12 CO(2-l). The centroid of each 
superbubble is shown with an asterisk. The position of the galactic center, which is the 
reference position of the maps, is shown as a cross. The range of channel velocities used to 
make each 0-th moment map is shown at the top right corner. The unit of intensity is Jy 
beam -1 km s -1 . The synthesized beam is shown at the bottom left corner, and the mosaiced 
primary beam is indicated by a dotted line at the 40% of its peak. The diamond in the 
SB1 map shows the position of the compact radio source called 5.17—45.4 in Ulvestad & 
Antonucci (1997). 
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Fig. 6. — Position-velocity diagrams across the superbubble SB1 in the direction parallel to 
the line of nodes (P.A. disk = 231°), (left) 12 CO(2-l), (right) 13 CO(2-l). The intensity is the 
excess brightness temperature in Kelvin, and is corrected for the primary beam attenuation. 
The rms noise in the data cube is 0.26 K for 12 CO and 0.10 K for 13 CO. Position offset is 
measured from the minor axis, while velocity is measured from the systemic velocity of 240 
km s _1 . The dotted ellipse shows the approximate location of the expanding bubble SB1 in 
each diagram. The black rectangle at the bottom left corner shows the spacial and velocity 
resolutions. The lowest contour is at 4 a for 12 CO and 3 a for 13 CO. The arrow in the left 
panel points the trough feature. 
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Fig. 7. — Position-velocity diagrams across the superbubble SB2 in the direction parallel to 
the line of nodes (P.A. disk = 231°), (left) 12 CO(2-l), (right) 13 CO(2-l). The intensity is the 
excess brightness temperature in Kelvin. Position offset is measured from the minor axis, 
while velocity is measured from the systemic velocity of 240 km s" 1 . The data are corrected 
for the primary beam attenuation. The dotted ellipse shows the approximate boundary of 
the expanding bubble. 
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Fig. 8. — Position-velocity diagram of the 13 CO(2-l)-to- 12 CO(2-l) ratio of excess brightness 
temperature. The p-v diagram is across the superbubble SB1 in the direction parallel to the 
line of nodes (P.A. disk = 231°). The ratio was calculated by using maps corrected for the 
primary beam attenuation and convolved to the same resolution. The ratio was calculated 
only where both lines are detected above 2.5 a. The black rectangle at the bottom left shows 
the spacial and velocity resolutions. The ratio approximately corresponds to the optical 
depth of the 13 CO(2-l) line. Note the low opacity gas in and around the superbubble, which 
is marked by an ellipse in a dotted line. 



